Eddy diffusivity in convective hydromagnetic systems 



M. Baptista^ S. M. A. Gama^ V. A. Zheligovsky^-^ 

^ Centre de Matematica da Universidade do Porto 
and Departamento de Matematica Aplicada 
Faculdade de Ciencias da Universidade do Porto 
Rua do Campo Alegre, 687, 4169-007 Porto, Portugal 
e-mail: mbaptist@fc.up.pt 

^ International Institute of Earthquake Prediction Theory and Mathematical Geophysics 
79 bldg.2, Warshavskoe ave., 117556 Moscow, Russian Federation 

^ Institute of Mechanics, Lomonosov Moscow State University 
1, Michurinsky ave., 119899 Moscow, Russian Federation 

^ Observatoire de la Cote d'Azur, CNRS 

U.M.R. 6529, BP 4229, 06304 Nice Cedex 4, France 



The date of receipt and acceptance will be inserted by the editor 



Abstract An eigenvalue equation, for linear instability modes involving large scales in a convective hydromag- 
netic system, is derived in the framework of multiscale analysis. We consider a horizontal layer with electrically 
conducting boundaries, kept at fixed temperatures and with free surface boundary conditions for the velocity 
field; periodicity in horizontal directions is assumed. The steady states must be stable to short (fast) scale 
perturbations and possess symmetry about the vertical axis, allowing instabilities involving large (slow) scales 
to develop. We expand the modes and their growth rates in power series in the scale separation parameter 
and obtain a hierarchy of equations, which are solved numerically. Second order solvability condition yields a 
closed equation for the leading terms of the asymptotic expansions and respective growth rate, whose origin 
is in the (combined) eddy diffusivity phenomenon. For about 10% of randomly generated steady convective 
hydromagnetic regimes, negative eddy diffusivity is found. 



1 Introduction 



According to the present-day paradigm, magnetic fields of most astrophysical objects - the Earth and the outer 
planets of the Solar system having a molten metal fiuid core [1], the Sun |2l and other stars, and even entire 
galaxies [3] - owe their existence to convective hydromagnetic processes y,l5|,|6|. Convection in the presence of a 
magnetic field obeys a familiar set of equations: the Navier-Stokes equation with Lorentz and Archimedes forces 
for the fiow, the magnetic induction equation for the magnetic field, and the heat equation for temperature. 
Virtually no analytic solutions of this system of equations are known except for some significantly reduced cases 
01 , under the assumption that certain symmetries are present, or for specific initial conditions. The referred set 
of equations may be used to simulate the evolution of astrophysical convective hydromagnetic systems. This 
approach was followed in @,&E3|5 where magneto-convection in an idealised plane layer was considered, and in 



nl . 12, 13, 14, 15, H, 17,ll8||, where the outer core of the Earth was modelled by three-dimensional equations 
of hydromagnetic convection in a spherical layer and, as a result, the predominant dipole morphology of the 
Earth's magnetic field was correctly reproduced in computations. 

However, in practice, accurate simulations for geo- and astrophysical real parameter values are close to im- 
possible, because the limited power of available computers prohibits computations with the adequate spatial 
(insufficient men iory) a nd temporal (insufficient CPU power) resolution. The simulations done by Glatzmaier 
and Roberts [111, [H HQ EH were performed for parameter values differing by several orders of magnitude 
from those characterising the outer liquid core of the Earth. Nevertheless, the agreement between these simu- 
lations and the geodynamo is surprising ^ISJ. The existence of sharp contrast spatial structures (e.g., Ekman 



boundary layer, which emerges in convective flows in a rotating Im^er with no-shp boundary conditions, and the 
instabihty of which may be a source of dynamo [ll, HO, 21, 22, 2^) and the prominent role which turbulence 



plays in the generation of magnetic flelds, indicate that such a high resolution is indeed necessary in simulations. 
Core-mantle coupling, which is believed to cause decade variations of the length of day, is another geophysical 
phenomenon involving small scales (topographic features at the boundary are unlikely to exceed 5 km in am- 
plitude; see [1]). Validity of numerical techniques for smoothening the resolution cut-off, such as employment of 
hyperviscosity, is questionable [24, 25, 26, 27]. 

A signiflcant uncertainty in rheology relations (28l| and parameter values [29] in the convective hydromagnetic 
equations (for instance, estimates of thermal diffusivity for the Earth core differ in several orders of magnitude 
jlj), makes it desirable to investigate typical regimes of behaviour of the solutions by varying the parameters in 
certain ranges and flnding locations, in the parameter space, of bifurcations marking drastic changes in behaviour 
(3Q| . a purely numerical approach to the implementation of this task would require a further signiflcant increase 
of the amount of computations. Consequently, a semi-analytic approach, employing asymptotic relations, is 
unavoidable as an alternative to both purely analytic and numerical approaches. Here we employ it to consider 
the problem of linear stability of three-dimensional convective hydromagnetic steady states in a layer. 

The characteristic spatial scale of the perturbed steady state is supposed to be much larger than that of the 
steady state. The ratio of the spatial scale of the flow (fast spatial variable) to the large scale of perturbation 
(slow variable), £, is a small parameter. (By small- and large-scale vector flelds we refer to flelds involving 
spatial scales of the order of the width of the layer and much larger scales^respectively.) Applying methods of 
the general theory of homogenisation of differential operators 31, 32, H, ll], we expand perturbation modes 



and their growth rates in asymptotic series in the parameter e and obtain a homogenised operator in slow 
variables, acting on mean flelds. Eigenvalues of this operator control stability to large-scale perturbations. The 
advantage of this approach stems from opening the possibility to disentangle the large and small scales and 
fully resolve small scales by solving the so-called auxiliary problems. 

Generically, the multiscale analysis reveals the presence of a-effect (see (H, 3^, 33|). The homogenised linear 
operator is then the flrst-order differential operator. Consequently, the system is generically unstable, since the 
spectrum of the operator is symmetric about the origin (if a mode W(x) is associated with an eigenvalue A, 
then W(— x) is a mode associated with the eigenvalue — A). In convective hydromagnetic systems which possess 
symmetry about a vertical axis or parity-invariance, a-effect is not present and the homogenised equations 
involve a second-order partial differential operator, whose eigenfunctions are Fourier harmonics. Its eigenvalues 
may be positive, implying instability. This phenomenon is referred to as negative (combined) eddy diffusivity 
[38|. Instability of this kind is weak: in the presence of a-effect the growth rate of the dominant perturbations 
is 0{e)^ whereas it is 0{e^) when a-effect is absent. Evaluation of eddy tensors emerging in the homogenised 
equations requires solution of auxiliary problems, which are linear elliptic partial differential equations in fast 
variables. With just a single characteristic spatial scale involved, they are not too demanding numerically. 

Multiscale asymptotic analysis was successfully applied to various problems o f hy drodynamics and magnetohy- 
drodynamics. The effect of negative eddy viscosity arises in two-dimensional [33, 4^, 41 1 and three-dimensional 
(42I, [il, Q hydrodynamic systems, if the flow is parity-invariant or if it is a Beltrami fleld (in j42| large scale along 
only one direction was assumed). Eddy diffusivity can be complex 45, 3|. In generic hydrodynamic systems, 
which do not possess the properties mentioned above, similar expansions indicate the presence of the so-called 
AKA-effect (i.e. anisotropic kinetic a-effect) jsl, [s^. In passive scalar transport systems, eddy diffusivity can 
only enhance molecular diffusivity (i^, [i^l . 

In the kinematic dynamo problem (concerning magnetic fleld generation, when the feedback influence of magnetic 
fleld on the flow via the Lorentz force is neglected), multiscale expansions were apparently flrst introduced in 
(iil and [42,] (where scale separation was related to fast rotation of the layer of conducting fluid). Similar 
asymptotic expansions in the kinematic problem (for flows, the amplitude of which may depend on the scale 
ratio) predict occurrence of a-effect [H^, 51, l^]. Generation of large-scale magnetic fleld by the negative magnetic 



eddy diffusivity mechanism is possible for parity-invariant steady [52|, |53|, iM] (in [M], large scale along only one 
direction was assumed) or time-periodic flows [55|, and by convective Bisshopp cell patterns [H^, symmetric 
about the vertical axis. Combined eddy diffusivity tensors for large-scale perturbations of both the flow and 
magnetic fleld constituting a parity-invariant three-dimensional MHD steady state were derived in [57]. 

In the papers cited above two different scales were present in the system. Multiscale expansions with three 
spatial scales were employed in [5^,[5§| to study the small-angle instability (g^I in convection in a rotating layer. 

Evolution of a mean hydrodynamic large-scale perturbation in the weakly nonlinear regime was considered in 
the absence of magnetic fleld for two-dimensional parity-invariant space-periodic flows |4ll. [6l| . and for three- 
dimensional MHD systems j62|. In j62| it is not required that the MHD state, nonlinear stability of which is 
examined, is either space periodic or steady; equations for the mean flow and magnetic fleld are generalisations 



of the Navier-Stokes and magnetic induction equation with an anisotropic (in general) combined eddy diffusivity 
tensor and quadratic eddy advection. 

In section O we present the equations for thermal convection in the presence of magnetic field and discuss 
boundary conditions and symmetries. In section [3l the multiscale formalism is applied. In section ISTH an eigen- 
function of the linearisation of a convective hydromagnetic system and its associated eigenvalue are expanded 
in a power series of the scale separation parameter and a hierarchy of equations is derived. In section 13. 2[ the 
solvability condition, which plays an important role in solution of equations of this hierarchy and in derivation 
of an equation for the mean part of perturbation in the leading order, is discussed. In sections [331 and [X4| the 
first two (order and order 1) equations in the hierarchy are expressed as a linear combination of the so-called 
auxiliary problems. In section 13. 5[ we consider the solvability condition for equations at order 2 and thereby 
derive the eigenvalue equation for the mean part of the leading terms in the expansions of the instability modes 
and their growth rates. At this stage emerges the homogenised combined eddy diffusivity operator acting on 
mean fields. In section [H we briefiy describe the numerical procedure for solving the auxiliary problems and 
present a set of basic fields which lead to large-scale instability for appropriate physical parameters (namely 
molecular diffusivities) . Finally (section [5]), we comment on possible extensions and limitations of the application 
of multiscale techniques employed here to study the instability of convective fiows in the presence of magnetic 
field. 



2 Equations of thermal convection in the presence of magnetic field 

2. 1 Time evolution of a convective hydromagnetic system 

Magnetic field generation by thermal convection is governed by the Navier-Stokes equation, the magnetic in- 
duction equation and heat transfer equation j63| : 

StV = V X (S X V) - + vd'^Y 

-H X (S X H) - a(T - To)G + F, 

a- V = 0, 

StH = 5 X (V X H) + T^S^H + R, (1) 
S-H = 0, 

dtT = -(V • d)T + kd^T + I |o> X H| V 5, 

where V = (^1,^2,"^), H = (i^i,i^2,^3) and T, depending on position in space x = {xi^X2^x^) and time 
t, are the velocity field, the magnetic field and the temperature, respectively. We use the notation dt = d/dt, 
di = d/dxi and d = Y^^=i ^idi, where e^ is the ith canonical vector. The term involving G (gravity, G = —ges 
for a horizontal layer) is the buoyancy force due to temperature variation and H x (9 x H) is the Lorentz 
force. F represents any other body forces acting on the fiuid, R is due to imposed external currents or magnetic 
fields, and S describes the distribution of external heat sources, u is the kinematic viscosity, r] the magnetic 
diffusivity, and a, k and a are parameters related to thermal expansion, thermal conductivity and electrical 
conductivity, respectively. Solenoidality of magnetic field follows from the Maxwell equations. Flows are deemed 
incompressible in line with the Boussinesq approximation. Henceforth, the system of equations ([T]) will be 
referred to as CHM (convective hydromagnetic). 

We consider the CHM equations in the spatial domain V = [0, Li] x [0, L2] x [0, L3], assuming periodicity in xi 
and X2 directions, and considering a finite layer in the xs direction. The boundary conditions at the surface of 
the layer are 

• for the velocity field: 

V3U3=0,L3 = 0? 

Ss^l 1^3=0, L3 = %V2|a;3=0,L3 = 0; 

• for the magnetic field: 

^3U3=0,L3 = 

93^^1^3=0, L3 = 1x3=0, L3 = 0; 



• for the temperature field: 



It is convenient to introduce the variable 

= T-Ti -STxs, 
where ST = (T2 — Ti)/Ls, satisfying the uniform boundary conditions 

^U3=0,L3 = 0- 

A solution to the CHM system may possess symmetry about the vertical axis, provided the forcing terms F, R, 5 
possess the same symmetry. A vector field Q is called symmetric (about the vertical axis) if 

Ql{-Xi, -X2,Xs) = -Ql{xi,X2,Xs), 
Q2{-Xi, -X2, X3) = -Q2{X1,X2, X3), 
Qsi-Xi, -X2,X3) = (53(^1,^2,^3); 

and anti-symmetric if 

Ql{-Xi, -X2,X3) = Ql{xi,X2,X3), 
Q2{-Xi,-X2,X3) = Q2{XI,X2,X3), 
Q3{-Xi,-X2,X3) = -(53(^1,^2,^3). 

We call a scalar field / symmetric if 

f{-Xi,-X2,X3) = f{xi,X2,X3), 

and anti-symmetric if 

/(-Xi, -X2,X3) = -f{xi,X2,X3). 

Symmetries are essential to eliminate first order (alpha) effects. In [43] parity-invariance is used to this purpose, 
but, for a horizontal layer, symmetry about the vertical axis is more realistic. Furthermore, parity invariance is 
inconsistent with the basic equations ([T]) for <j 7^ 0. Under an appropriate forcing, any hydromagnetic convective 
system will possess steady states with these symmetries. 



2.2 Linearised CHM operator 



Let us consider a steady state solution, p, V, H and ^, of the CHM system ([T]) and a small perturbation, pe'^^ 
Ve^*, He'^^ and Oe^^ ^ of this steady state, where p, V, H and depend only on spatial variables. In what 
follows, we will call the spatial profiles of the perturbation fields, p, V, H and 6>, a perturbation. Replacing p, 
V, H, 6>, respectively, by p + pe^^ V + Ve'^^ H + He^^ 9 + Oe'^^ in ^ and neglecting second order terms in 
the perturbation, we obtain an eigenvalue problem for the perturbation: 

dp 

AW = AW 



v = o, 

S-H = 0. 

Here, the block notation introduced in ji^ is used: 



(2) 



W 



V 
H 



(3) 



is the (3 + 3 + l)-dimensional block column vector combining the 3 components of the fiow, the 3 components 
of the magnetic field and the temperature field. In what follows, (3 + 3 + l)-dimensional vectors of a similar 
structure, will be used. The operator A is obtained by linearisation of the CHM equations in the vicinity of the 
steady state V, H and 0] it can be represented as a block matrix (acting on (3 + 3 + l)-dimensional vectors 
of the structure similar to (O): 



(4) 



+ V X (S X •) - (S X V)x -H X (S X •) + (a X H)x -aG 
X (H X •) 7^52 + 5 X (V X •) 

_-(• • d)0 - STe3' cr(S X H) • (S X •) kd^-V-d 

Note that A preserves the symmetry of fields W, symmetric (or anti-symmetric) about the vertical axis. 

The complete formulation of the eigenvalue problem ([2]) involves specifying spatial periods of perturbations 
which can be any integer multiples of the periods Li and L2. If the smallest of the periodicity boxes is considered 
the system of equations ([2]) is referred to as the problem of linear stability to short-scale perturbations. 



3 Linearised large-scale CHM equation 



3.1 The two-scales expansion 



In this section we construct a homogenisation of the hnearised CHM operator. Eigenvalues of the homogenised 
operator control linear stability of the CHM steady state to perturbations with spatial periods large enough for 
the asymptotic behaviour to set in. 



Following the method applied in previous studies [52|, |53|, l57|, l55| , we consider fast variables, x, representing the 
short scale dynamics, and slow variables, X = (Xi,X2) = (5x1,6X2), representing the large scale dynamics. (In 
a layer of finite width only slow variables in horizontal directions are geometrically consistent.) The parameter e 
is the scale separation. The perturbations p, V, H and 9 are assumed to depend on both fast and slow variables. 
We denote Vi = d/dXi and V = Yld=i ^i^i- 

Now we perform differentiation in ([2]) applying the chain rule: di ^ di -\- eVi for z = 1, 2. This implies 

A = A(0)+eAW+e2A(2), 



where 



A(o) 



Ad) 



i/S^ + V X (5 X •) - (5 X V)x -H X (S X •) + (S X H)x -aG 



X (H X •) 



rjd'^ + 9 X (V X •) 
a{d X H) • (S X •) 





kd^-V-d 



2zyS • V + V X (V X •) -H X (V X •) 
- V X (H X •) 27^S • V + V X (V X •) 

cr(S X H) • (V X •) 2/c5-V-V-V_ 



A(2) = sv^ with H 



T] 

0k 



(5) 



(6) 



(7) 



(H is the molecular diffusivity tensor). Note that A*^^) and A*^^) preserve the symmetries of both symmetric and 
anti-symmetric fields, but A^) exchanges their symmetry. 

Next we expand W, p and A in power series of e: 

n 



i=0 



(8) 
(9) 



(10) 



i=0 



Each coefficient W^'^ and p^*^ in the expansions is a function of both x and X. Substituting these series in (0) 
and equating the terms in e", at each order n, we obtain a hierarchy of equations: 



order 0: 



order 1: 









a-v(°) = 0, 

5-B(o) =0, 



A(o)w(i) = -A(1)W(0) + AoW(i) + AiW(o) 









(11) 

(12) 
(13) 



(14) 

(15) 
(16) 



— order 2: 



a.v(2) = -V-VW, 
5.b(2) = -V-B«, 



ap(2) + ^p(l) 







order 



i=0 









5.b(") = -v-b("-i). 



(17) 

(18) 
(19) 



(20) 

(21) 
(22) 



Let (•) = {LiL2Ls)~^ •dxidx2dx^ denote the mean (over the fast variables in V) and {•} = • — (•) denote 
the fluctuating part of a vector or scalar fleld here denoted indistinguishably by •. The average (•) is the large 
scale component of the respective fleld. It is possible to solve recursively all equations in the hierarchy, flnding 
all terms of the expansions (|8))- (p!Q)) . Only the equations up to order 2 are required to derive a homogenised 
eigenvalue equation for the mean parts, (W^^^) and (p*^^^), of the leading terms. It emerges as the solvability 
condition for the equation in the fast variables at order 2. 



3.2 Solvability Conditions 



Let P be the projection into the subspace of (3 
components are solenoidal: 



3 + l)-dimensional flelds, in which the 3-dimensional vector 













p 





















where d'^Q^^ = d • and d^Q^^ = d • Q^. In what follows, we will solve equations of the form 

PA(o)f = p^, 

where vector components of f are required to be solenoidal; this equation is thus equivalent to 

PA(o)p/ = Pg, (23) 

By the Fredholm alternative [g^], a solution of ([23)) exists if and only if is orthogonal to the kernel of 
(PA'^^^P)*, where * denotes the adjoint operator. In other words, the solvability condition for ([23)) is (P^, c) = 0, 
where c is any vector in ker(PA'^^^*P) and (•,•) denotes the C'^ inner product. As usual, the adjoint operator 
can be derived performing integration by parts in the identity (A*^^^* Wi, W2) = (Wi, A*^^^ W2). In the 
present case 



i/S^ - S X (V X •) + (S X V)x H X (S X •) -esST + Od 

S X (H X •) - (S X H)x rjd'^ -V x{dx •) x (S x H) - a{d x H) x {^^) 

-aG ■ kd^^V -d 



Boundary conditions for the vector flelds in the domain of A^^^ can be found from the condition that boundary 
surface integrals, emerging in integration by parts in the scalar product (Wi, A*^^^ W2), vanish. It can be verifled 
that the boundary conditions that we assume hold for vector flelds in the domain of the adjoint operator as 
well. 

Our construction relies on the existence of vector flelds in ker(PA^^^P) with non- vanishing average horizontal 
components of the flow and magnetic fleld (see a detailed discussion in [56]). It can be easily established that the 
dimension of the subspace of ker(PA'^^^P), consisting of such vectors, is equal to the dimension of the subspace 



of ker(PA*^^^*P) consisting of vectors with non-zero horizontal space averages. If c = 



where and 



are any constant horizontal vectors, evidently PA'^^^^Pc = 0, since A'^^^ c is a gradient. Such constant 
vector fields satisfy the boundary conditions under consideration; therefore any such c belongs to ker(PA'^^^*P). 
Thus, ker(PA'^^^P) is at least four-dimensional. In what follows, we assume that the dimension is four, which 
is generically the case, and thus any non-zero vector from ker(PA'^^^P) has non-zero horizontal averages of the 
fiow and/or magnetic field components. Then the solvability condition for ([23)) consists of orthogonality of 
and to constant horizontal vectors, i.e. the horizontal components of (g^) and (g^) must vanish. 



3.3 Equations at order 



Decoupling of the large- and short-scale behaviour is evident in the form of solutions of the equations emerging 
at orders and 1 in the hierarchy. Since A*^^^ is an operator in the fast variables, short-scale variation of these 
solutions is exclusively due to their multiplicative dependence on solutions of ([TT|) - ([T3)) and ([T4)) - ([T6)) - the 
so-called auxiliary problems, which are partial differential equations in the fast variables. 

Performing integration and using the boundary conditions, we find that the horizontal components of the fiow 
and magnetic field components of (A'^^^W'^^^) vanish. Since we seek solutions where (V*^^^) and (H*^^^) do not 
vanish simultaneously, averaging of (pT|) implies Aq = 0. The problem to solve becomes 









By linearity. 



The fields 



and are linearly independent solutions of the problems 



(24) 



(25) 



a-sv = o, 
a-s" = 0, 









(26) 



(27) 



i.e. Si are linearly independent fields in ker(PA'^^^P). The gradient part of ([26)) satisfies 

s^^f = s- (a(o)s,)''. 

Note that {p^^^) cannot be determined at this order, since derivatives in the fast variables in the r.h.s. of ([24)) 
eliminate any averages in p^^\ Si and depend only on the fast variables. The coefficients of the linear 
combinations f ([25)) ). a^, depend only on the slow variables. 

The four possible averages (constant mean fields) in ker(PA*^^^P), selected by the boundary conditions, are 





ei 




62 


(Si) = 












_ _ 




_ _ 




" " 




" " 


(S3) = 


ei 








_ _ 




_ _ 



Hence, the four auxiliary problems at order are 

A(°){Si} = 

A(°){S2} = 



5(5f - Vi) + aiV 

a(5f - V2) + a2V 
d20 



(28) 



(29) 



AC) {S3} 



A(0){S4} 



-disi 



'd{Sl 



■Hi)- 
-dx\ 


-H2)- 




diR' 



(30) 



(31) 



with Si given by (|27l) . 

Solenoidal parts of the right hand sides of (|28)) - (|3T|l are anti-symmetric and, since A^^^ preserves the symmetry 
of fields, Sj are anti-symmetric. Evidently, the order auxiliary problems are of the form of ()23)) . and their 
solvability follows immediately from the periodicity of the CHM steady state V, H, 6 in horizontal directions. 



3.4 Equations at order 1 
Averaging of (IHl) yields 



A(i)w(o)) + Ai{W(o)) 



V(p(o)) 





(32) 



As in order 0, horizontal parts of the flow and magnetic fleld components of (A*^^^ W*^^^) are zero. The same holds 
for (A'^^^W'^^^) (this can be shown, integrating directly the terms of the form of spatial derivatives and using 
the boundary conditions, and exploiting the symmetry of the perturbed CHM steady state and anti-symmetry 
in the fast variables of V^^^ and H^^^ when considering mean horizontal parts of the terms of the form of 
V X (V X V(o))). (|32]) becomes 

Ai(W(o)) 

Thus, if Ai ^ 0, then (H*^^)) = and (V'^^)) = -V {p'^^^) / \i. However, averaging of (dH) over the fast variables 
yields V-(V'^^^) = 0, with (V*^^^) and V{p^^^) belonging thereby to orthogonal subspaces. This implies (V^^^) = 0, 
which contradicts the original assumption that (V^^^) and (H*^^^) do not vanish simultaneously. Therefore, 
Ai = 0, {p^^^) = and (dH) reduces to 







0" 



















'vp^^y 









+ 














(33) 



From 



where 



and 



we flnd 









i=l i=l 









2vdj 



He^ • 



-e^H ■• + Hj 












(34) 



The problem at this order reduces to 
Therefore, by linearity. 



4 2 

i=i j=i 









4 2 4 

1=1 j=i 



(35) 



2 = 1 



i=l j=l i=l 



(36) 



where Fy 



and rP- are mean-free linearly independent solutions of the auxiliary problem at order 1 : 



A(0)ry = My 



d-ry 









(37) 



Taking the divergence of the velocity component, we obtain a Poisson equation for Ff-: 

V 



a2/;p=a.(A(°)ry-My) . 



(38) 



The average of p^^^ cannot be determined at this order, since derivatives in the fast variables in the r.h.s. of 
(|33)) eliminate it. In (|35|) and (|36)) . hi depend only on the slow variables, and the fields Tij and Ff- only on the 
fast ones. It is convenient to solve (|37|) in the subspace of (3 + 3 + l)-dimensional vector fields, where vector 
components are solenoidal. Consider the substitution 







The conditions 



imply d • T'Yj — ^ ' ^'Tj = 0- At order 1 we have thus to solve eight equations: 



/H 



PA(o)r',-, = p Mi, - a(o) 







(39) 



for i = 1, . . . , 4 and j = 1,2. Solvability of (|39)) can be easily verified by symmetry arguments, since Bj changes 
the symmetry of fields. 



3. 5 The mean-field equations for the CHM instability mode 



At order 2 the solvability condition is not trivially satisfied and yields equations for the large-scale mean 
components of the instability mode. We consider orthogonality of the r.h.s. of (p!7|) to ker (PA*^^^*?), i.e. 



(Ci, A2W(°) - A(2)w(o) - A(i)w(i) 

VC( e ker(PA(0)*P). From 









0, 



(40) 



4 2 

i=l j,k=l 



Since Ci are constant, (|4Q)) is equivalent to 



4 2 



(41) 



From the system of equations (|4T]) . we find a^; then W*^^^ is obtained from ([25|) . Thus, we have derived a closed 
set of equations for the leading terms in the expansions (|8))- (p!Q)) of eigenmodes and their growth rates. The 
leading term in the eigenvalue expansion is A2, i.e. A = 0{e'^). This growth rate determines the characteristic 
slow time scale of the large-scale dynamics: T = e^t. 



Since only horizontal components of and can be nonzero for constant vectors G ker(PA'^^^*P) (see 
section [321) , we can choose Ci = (Si). This implies (S^, C^) = Su (here Su is the Kronecker symbol). Then (|4T]) 
takes the form 



i=l i=l j,k=l 

+(cr,v(p«)) = o. 



(42) 



This is an eigenvalue problem for the second order partial differential operator with constant coefficients, which 
is called combined eddy diffusivity operator. It admits Fourier harmonics as eigenfunctions: 



a„(X) = a„(q)e*'i-^, (p' 



p(q)e 



iqX 



(43) 



where q = (<7i, ^2) is an arbitrary unit wavevector and n = 1, 4. Upon substitution we find that the coefficients 
dn satisfy 



(44) 
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where E is the 4x4 matrix 



E/i = ^ qkqjiCi.BkTij). 



(45) 



Averaging of ([T5)) and ([T6)) yields V • (V^^^) = and V • (H*^^^) = 0. By virtue of these solenoidality conditions 
and gni), 

(ai,a2) = d[{q2,-qi), 
{ds.d^} = a2((72, -qi)' 

Substituting these expressions into (|44|) and scalar multiplying it by (^2, — ^1; ^2, — ^1), we reduce (|44|) to an 
equivalent 2x2 eigenvalue problem: 



(A2 + iy)d[ 
(A2 + 7/)a2 



E' 



= 0, 



(46) 



where 



E^i = Eng^ _ (E12 + E2i)gig2 + E22gi, 

^12 = Eia^l — (Ei4 + E23)^1^2 + E24^i, 
^21 = Eai^l — (E32 + ^4i)QiQ2 + ^42Qij 
^22 = Eaa^l — (E34 + ^43)QiQ2 + E44g^. 
Noting that q = (cos 6>, sin ^), 6 G [0,27r], we obtain 



A±(^) = -- l±^/l-^ 



4c 
6 



with 6 = ly + 7] + E'li + E'22 and c= vrj-\- vE'22 + ??-£'n + E'iiE'22 — ■£'i2-E'2i- The maximum and minimum growth 
rates, 



A5 



\ mm 
^2 



max max 

6'G[0,27r] 



are admitted for 6>'s denoted by 6>^^^ and 6^' 



respectively. 



mm mm 

6'G[0,27r] 



(47) 
(48) 



4 Numerical Results 



The auxiliary problems were solved numerically using pseudo-spectral methods to evaluate the action of the 
operators A*^^^ ([5|) and Bj (|34|) on the fields. In the finite direction of the layer, the usual plane wave basis was 
replaced by a half period sine or cosine basis, agreeing with the boundary conditions: 

for a scalar function satisfying Dirichlet-kind boundary conditions, and 

J nky , rik^ 

for a scalar function satisfying Neumann-kind boundary conditions, with kx = 27rn/c^/Li, ky = 27rn/c^/L2, 
kz = TTTik^/Ls and n/c^, n/c^, G Z. For each auxiliary problem, a linear system of equations in the Fourier 
space was obtained and solved numerically by the conjugate gradients method [65]. 

Asymptotic expansions for large molecular diffusivities, as well as comparison with previous calculations for plan 
form velocity fields [56], were used to validate the code. As previously stated, the basic steady state must be 
stable to short-scale perturbations, i.e. the dominant eigenvalue (Xshort) of the operator A must have a negative 
real part. The dominant eigenvalue can be evaluated using the method used in [66] for perturbations in each of 
the two symmetry subspaces. 

We want to model magnetic instabilities in turbulent convective fiows. For the reasons exposed in the intro- 
duction, simulations of fully turbulent regimes are very resource expensive. Within the scope of our approach, 
steady states can be randomly generated with decaying energy spectrum. Such states satisfy the basic equations 
for the appropriate source terms. Usually only a finite number of Fourier harmonics {kmin < k < kmax) is 
generated, the remaining being set to 0. Applying the appropriate linear transformations, we make sure that 
the generated CHM steady states are solenoidal and possess the required symmetry. The coefficients are then 
normalised to obtain the desired energy spectrum and the norm (r.m.s.) of each field is set to 1. Algebraic 
{E{k) ^ k~^) or exponential {E{k) ~ exp(— ^/c)) spectra were used in (Hsl. IHtI. . 

Zheligovsky et. al. found that fiows with exponentially decaying spectra are statistically better dynamos. 
However, in fully developed turbulence, the energy spectrum in the inertial range is known to be algebraic. In 
this work, all fields have been normalised to have decaying algebraic energy spectra, with ^ = 4, for the Fourier 
modes with < A: < 7. Simulations have been carried out for the periodicity box of size 27r x 27r x tt, with the 
resolution of 32 x 32 x 16 Fourier harmonics. An ensemble of 1000 instances of CHM states has been generated. 
It turns out that 110 out of 1000 generated fiows exhibit negative combined eddy diffusivity (see Fig. [T]). The 
values ly = jj. = K, = O.b were chosen that large so that to make sure that the randomly generated CHM states 
were stable to short-scale perturbations. We have directly checked, by computation of the decay rates of the 
dominant short-scale modes, that 30 of the generated CHM states from our ensemble are indeed stable; for 3 of 
them eddy diffusivity is negative. No instances of CHM states unstable to short-scale perturbations were found 
for these values of diffusivities. 
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Edd^ Diff jsiuity 

Figure 1 Statistics of eddy diffusivity. 



In figure [21 the fields corresponding to one of the generated steady states are presented. The short-scale growth 
rates are Xghort = —0.5662, for symmetric perturbations, and Xshort = —0.05175, for antisymmetric perturba- 
tions. The maximum and minimum growth rates of large-scale perturbations are X^"^^ == 1.426, for 0'^^^ = 3.392, 
and A2^*^ = —1.165, for 6^'^^ = 3.593, respectively. The maximum growth rate is positive (which corresponds 
to a negative eigenvalue of the eddy diffusivity tensor), i.e. a large-scale instability is present. All the auxiliary 
problems show a decaying energy spectrum (see figures [MH] below) and the expected symmetries can be observed. 

Spatial short-scale structure of the large-scale eigenmode is defined by the leading term W*^^^ in the expansion 
dl]), and therefore by the fields (see ([25]) ). Magnetic field in all of them (see figures [MB]) has the form of 
distorted "retrograde columns" [67] and concentrates near the horizontal boundaries of the layer. This kind 
of structure is reproduced in their linear combination W*^^^ for the most unstable mode (see figure [7|). Such 
behaviour was originally noticed in the non-linear evolutionary simulations [67] of magnetic field generation 
by rotating thermal convection, and it can be attributed to the boundary conditions, namely, ideal electric 
conductivity of the boundaries. It is interesting that this feature is quite robust - from the formal point of view 
our problem is evidently quite different from the one considered in [67] . 
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Steady state temperature Steady state energy spectra 

Figure 2 Basic (steady state) fields and corresponding energy spectra. Symmetry about the z-axis is observed in all of 
the displayed fields. 
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S4 : temperature S4 : energy spectra 

Figure 6 S4 fields and corresponding energy spectra. Anti-symmetry about the 2;-axis is observed in all of the displayed 
fields. 





5 Concluding Remarks 



We have derived an eigenvalue equation for large-scale perturbation modes of a CHM steady state. On average 
(over small spatial scales) the modes are, in the leading order, simple harmonic waves. Their growth rates are 



controlled by the combined diffusivity tensor, involving molecular kinematic viscosity and magnetic diffusivity 
and an additional tensor - the so-called combined eddy (turbulent) diffusivity correction, which is anisotropic 
(the entries of the matrix depend on the direction of the wave vector q), and which intermixes the influence 
of the flow and magnetic fleld. Originally this mutual influence is due to advection (the influence of the flow on 
magnetic fleld) and the action of the Lorentz force (the influence of the magnetic fleld on the flow), but it has 
now a different algebraic form - in particular, unlike in these basic physical laws, it is, on average, linear. 

We have found that about 10% of randomly generated steady CHM regimes, that are stable to short-scale 
perturbations, exhibit negative eddy diffusivity; such steady states are unstable to large-scale perturbations. 
However, the growth rate of the perturbation is quadratic in the scale ratio e, i.e. it is small. Thus, this instability 
can be observed only if the considered CHM steady state is stable to short-scale perturbations, which would 
have larger {0{e^)) growth rates otherwise. Other competing linear instabilities may also persist. For instance, 
in thermal convection in a rotating layer with free boundaries (with no magnetic fleld present), steady rolls, 
steady square cells, standing and travelling waves near the onset were demonstrated [58*, ^ to be unstable 
to large-scale perturbations of a particular form (the small angle instability), with the growth rates scaling as 
0(6^) for the rolls and the cells, and 0{e) for the waves (here £ is the smallest scale ratio in the system); the 
rolls and square cells possess symmetry about the vertical axis. 

The restriction that the perturbed CHM state is steady can be lifted [62*|. Instead of averaging over fast spatial 
variables, averaging over the spatio-temporal domain of fast variables must then be performed. This allows, in 
particular, to carry out the stability analysis of time-periodic CHM states. However, in any case the perturbed 
states are required to be symmetric (for instance, symmetry about the vertical axis, as considered here). Parity 
invariance is another type of symmetry consistent with the CHM equations, if the Joule term is neglected. Like 
the symmetry about the vertical axis, it guarantees that no a-effect emerges and it allows to construct a second 
order combined eddy diffusivity operator by the same method of homogenisation. 

Different multiscale expansions are obtained for different sets of conditions imposed at the horizontal boundaries 
of the layer. Often considered alternative boundary conditions include the no-slip condition for the flow, an 
isolator outside the fluid layer condition for the magnetic fleld, and heat insulating boundaries (zero heat flux 
condition) for temperature. They can be of interest, for instance, in geophysical applications (for which the no- 
slip condition at the outer kernel boundaries and the isolator condition at the outer boundary of the spherical 
layer are more appropriate than those considered here) . The method of homogenisation that we have used relies 
on the existence of constant vector flelds in ker(PA'^^^*P). In addition to the two constant horizontal vectors in 
the flow and magnetic fleld components considered here, another scalar constant, a flxed temperature, belongs 
to ker(PA*^^^*P), if Joule heating is neglected (a = 0) and the zero heat flux condition for temperature is 
assumed. Therefore, the procedure of homogenisation that we have used can be applied if, at least, one of the 
following conditions is imposed on the horizontal boundaries: free boundaries, or conducting boundaries, or 
no heat flux. Each quantity, satisfying boundary conditions from this list, increases dimension of the problem 
for large-scale mean-flelds obtained from the solvability condition for the equations emerging at order 2. The 
remaining ones (temperature in the particular case that we have considered here) are essentially short-scale and 
only affect, via solutions of auxiliary problems at orders and 1, the values of coefficients in the combined eddy 
diffusivity tensor for the large-scale quantities. The boundary conditions for which this approach is not directly 
applicable (i.e. the no-slip boundary condition for the flow, or the insulator condition for the magnetic fleld, 
or flxed temperature) can apparently still be treated by the homogenisation method that we have applied, but 
this requires considering boundary layers, increasing signiflcantly the complexity of the problem. 

A common feature of astrophysical convective systems, such as interiors of planets or stars, is rotation. A 
straightforward incorporation of the Coriolis force in the analysis is inconsistent with the homogenisation pro- 
cedure that we have used: averaging of the linearised Navier-Stokes equation (for the free boundary conditions) 
shows that constant average non-zero horizontal velocities give rise to a non-zero average Coriolis force, which 
can be balanced only by a constant gradient of pressure. This suggests an unbounded linear growth of pres- 
sure, which is not, however, unphysical: in a rotating system only pressure can offset the centrifugal force. The 
simplest way to overcome the resultant algebraic difficulties is to consider the vorticity equation, for which the 
methods for construction of the two-scale expansion are applicable without any modifications. 
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